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Abstract 

We investigate the autocorrelator of conductance peak heights for quantum 
dots in the Coulomb bfockade regime. Analytical and numerical results based 
on Random Matrix Theory are presented and compared to exact numerical 
calculations based on a simple dynamical model. We consider the case of 
preserved time-reversal symmetry, which is realized experimentally by varying 
the shape of the quantum dot in the absence of magnetic fields. Upon a proper 
rescaling, the correlator becomes independent of the details of the system 
and its form is solely determined by symmetry properties and the number of 
channels in the leads. The magnitude of the scaling parameter is estimated 
by a semiclassical approach. 

PACS numbers: 73.20.Dx, 05.45. +b, 72.20.My 
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I. INTRODUCTION 



Since the first experiments deahng with coherence phenomena in microstructures in 
the 80's and the theoretical activity that followed, many experiments have probed several 
different universal features of mesoscopic devices. One important example of such systems 
are semiconductor quantum dots the subject of our study. With the fast progress in 
nanotechnology, the experiments gained a lot in sophistication and were able to explore 
more subtle aspects of mesoscopic physics. One of the most interesting subjects in the field 
is related to the influence of chaos in the ballistic electronic transport p|,|^,§J|,|,^||,|| . 

Presently, we know that the statistical measures of conductance (distributions and cor- 
relators) in open systems are well described by universal curves obtained from Random 
Matrix Theory (RMT). However, the understanding of system-dependent properties (aver- 
age conductance, correlation scales of external parameters, etc.) demands a complementary 
approach, which extracts relevant elements of the underlying classical dynamics or average 
quantum properties of the system at hand. In particular, the decay width of the conductance 
correlators as a function of either an external magnetic field or of variations in the shape of 
the quantum dot can be understood in semiclassical terms. 

In this work we extend our recent analysis of the correlator of Coulomb blockade 
conductance peak heights for quantum dots in the chaotic regime. Here we concentrate 
on the case of preserved time-reversal symmetry and generalize our previous semiclassical 
analysis. We use the same notation as in our previous paper. 

The conductance G of mesoscopic devices, such as quantum dots, can be evaluated using 
the Landauer-Biittiker formula ||To|JTl[] . For quantum dots the latter becomes particularly 



simple and can be expressed in terms of the partial decay widths Vcv and partial decay 
amplitudes '-)cu 0: 




^cv = M , with = J— dsx*c{r)^y{r) , (1) 



where V'l/ is the eigenfunction of the resonance v inside the dot with the appropriate boundary 
condition and Xc is the wave function at energy for the channel c belonging to one of 
the leads. The integration is performed over the contact region between the lead and the 
quantum dot. Equation (|l]) does not include barrier penetration factors 0, which should not 
be important for the case of one open channel in each lead - the case we analyze here and the 
one which is presently more relevant for comparison with experiments. In this framework, 
the influence of finite temperatures caused by the rounding of the Fermi distribution can 
be done straightforwardly. In the single-level limit T^v <^ kT < A (A is the single-particle 
mean level spacing), the conductance peak height corresponding to an on resonance 
measurement is given by |12[ 



where Tii^^^y is the width of the resonance v corresponding to a decay into open channels in 
the L{R) lead. In other words, Vk^r^ = I]ce-L(H) ^cu- If the quantum dot has any asymmetry 
(or if there is a small amount of disorder) and the leads are placed far from each other on 
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the scale of the electron wavelength, and Tlu will not be identical and should, in fact, 
vary independently. 



II. PARAMETRIC CORRELATORS AND RANDOM MATRIX THEORY 

We are interested in universal properties of a system Hamiltonian H which depends on an 
external parameter X and whose underlying dynamics is chaotic. Without loss of generality, 
we shall illustrate our study by considering X as a measure of the shape deformation of the 
quantum dot boundaries (in Ref. the parameter X was identified with the magnetic flux 
threading the quantum dot). Furthermore, we assume that in the absence of a magnetic 
field, H can be modeled as a member of the Gaussian orthogonal ensemble (GOE) of random 



matrices [|T3|. Thus, the statistical properties of the system should be obtained by an 
ensemble averaging. We focus our study on two autocorrelation functions: the level velocity 
correlator Ct,(X) and the conductance peak correlator Cg{X). 



As originally found by Szafer, Altshuler, and Simons |ll^JT5[| the level velocity correlator 

'de^{X-X/2) de^{X + X/2)\ / de^{X)\^ 



CM ) - ^2 



dX dX I \ dX I 



(3) 



is a universal function for systems whose underlying classical dynamics is chaotic. Here (■ ■ ■) 
denotes an average over resonances v and over different values of X. Due to the ergodic 
nature of the systems under investigation, one usually conjectures that (■ ■ ■) is equivalent 
to an average over an ensemble of random Hamiltonians. The universality becomes explicit 
after applying the rescalings 



X = X^a(O) and c„(x) = C^{X)IC^{<S) . (4) 

It is very difficult to measure Cv{X) experimentally for quantum dots and, consequently, to 
obtain the scale C„(0). The main goal of this work is to address this difficulty by suggesting 
an alternative way of studying parametric correlations. We show in Section IV that one 
can estimate C„(0) by semiclassical arguments, once details of the confining geometry of the 
dot are known. An important statement implicit in Refs. [T^ and is that the quantity 
^Jcjfi) sets the scale for any averaged parametric function (/(X)). Therefore, we also study 
another (experimentally more accessible) function, the conductance peak height correlator 
Cg(X), which is defined as iJI 

C,iX) ^ (g^ - X/2) + X/2)) - {gu{X)f . (5) 

This correlator also becomes universal upon rescaling X according to Eq. (|) and Cg{x) = 
Cg{X)/Cg{0). For the GOE case an asymptotic expansion of Cg{x) for x <^ 1 can, in 
principle, be found analytically, in analogy to the GUE case shown in Ref. [Q. Here, however, 
we only present results based on numerical simulations, which nevertheless cover the entire 
X range. For that purpose, we have performed a series of exact diagonalizations of random 
matrices of the form H{X) = Hi cos{X) + H2 sin(X), with Hi and H2 denoting two 500 x 500 
matrices drawn from the GOE. This model and the details of the calculation are discussed 
in Ref. 0. The final results for c^(x) and Cg{x) are presented as full lines in Figs. |l| and |^, 
respectively. 



3 



III. DYNAMICAL MODEL 



The aim of this section is to compare the statistical resuhs of the previous section with 
exact numerical diagonalizations of a dynamical model. The essential characteristic of a 
dynamical model for this type of study is a fair resemblance with the actual experimental 
conditions, combined with its adequacy to numerical computations. For this reason we chose 
the (two-dimensional) conformal billiard [1T7|,|5|,|9[| . Using complex coordinates, the shape of 
the billiard in the w plane is given by |2| = 1 in the area-preserving conformal mapping 
w{z) = [z + hz"^ + ce'^^z^)/\/l + 26^ + 3c^, where h, c, and 5 are real parameters chosen in 
such manner that > for \z\ < 1. In particular, we chose the values h = c = 0.2 

and obtained different boundary shapes by varying 5. Hence, X is now to be identified 
with variations in 5. The classical dynamics of a particle bouncing inside billiards defined 
by these parameter values is fully chaotic 0. The quantum dynamics is determined by 
choosing Neumann (zero current) boundary conditions, which yield non-trivial behavior at 
the boundary, as needed in Eq. (p. The eigenstates ipu thus obtained correspond to the 
resonant wave function appearing in Eq. (||) |]5|. 

In practical calculations, one uses as a truncated basis the lowest (in our case 1000) 
eigenstates of the circular billiard (6 = c = 0), which are known analytically. To solve the 
Schrodinger equation in this basis one has to calculate several thousand matrix elements 
of the Jacobian J = \w'{z)\'^. However, using polar coordinate separation, only matrix 
elements of powers of the radial coordinate need to be calculated numerically. The angular 
part is done analytically and changes in shape are trivial, since b, c, and 6 act as prefactors 
to the matrix elements. We have calculated the spectrum for 0.41 < 6 < 2.73 in steps of 
0.04, whereby the two values 5 = 0, vr which lead to spatial symmetry were avoided. Due 
to the truncation of the basis, only the lowest 250 of the calculated 1000 eigenstates were 
accurate enough to be used in the analysis. We also discarded the lowest 150 eigenstates 
because of their markedly nonuniversal behavior. The spectrum was then unfolded using 
the Weyl formula. We emphasize that it is not a trivial task to increase the number of 
converged states. No symmetry reduction of the resulting eigenvalue problem is possible for 
the asymmetric conformal billiard. 

The scaling parameter C„(0) was calculated using Eq. (^ and averaging over groups of 20 
levels and 12 shapes. For magnetic flux variations we had found that (7^,(0) is proportional 
to e^/^ 0. Here, for shape variations, we find the following fit: 

CMe,5 = as e^/\ (6) 

with 6 = 0.0003 + 0.00145. A physical explanation for this energy dependence will be given 
in the next Section. We do not have, however, a full undestanding of the 5-dependence of 
as yet. Note that the numerical coefficients in as are specific to a billiard with an area equal 
to vr^. 

We found that spectral fiuctuations were rather large and for that reason it seems more 
difficult to use shape than magnetic flux as the external parameter to obtain good statistics, 
at least in numerical simulations. Because of the ^-dependence of Cy{0), the rescaling of X 
has to be performed with care. The linear relation x = ^JCv{0)X holds only locally. We 
therefore used the generalized relation 
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x{e.X)= J^"" ^CMe,s d6. (7) 

In Figs. ID and|^ we show the numerical results for Cy{x) and Cg{x), respectively. They both 
compare rather well with the curves obtained from the GOE matrix diagonalizations. 



IV. SEMICLASSICAL ESTIMATE OF 0^(0) 



In a recent paper [0, Berry and Keating obtained an approximate expression for the 



level velocity correlator [Eq. (y)] by combining derivatives of the cumulative level density 
expressed semiclassically with the Gutzwiller trace formula. Their calculation applied for the 
case where X is an Aharonov-Bohm flux line. In what follows we generalize the derivation 



presented in Ref. [|T6[ for an arbitrary parametric variation of the Hamiltonian. Furthermore, 
here it will become more transparent why the result is universal for any chaotic system, 
depending only on whether the time-reversal symmetry is preserved or not. We do not 
consider the case of broken rotation symmetry (symplectic ensemble). 



The derivation presented in Ref. |]16[ approximates the velocity level correlator as the 



following function of an external parameter X: 

where the average is taken over a window of width SE ^ A around E. One should view this 
approximation with a certain caution. It is in fact exact for X = once the regularization 
parameter t] is gauged properly. For large values of X, Eq. contains the leading order 
term of Cy{X). For any other (intermediate) values of X, F^ just interpolates smoothly 
between the correct asymptotic values. The inaccuracy over intermediate values of X does 
not rely a priori on the approximation used in the evaluation of F^, but has a trivial origin: 
One cannot express C„(X) exactly by a 2-point, fixed energy correlator, as in Eq. (H). 

Equation (|^) can be evaluated semiclassically by expressing the cumulative level density 
A^^ by the trace formula Njj{X; E) = J2fi exp(?S'^/^ — r]T^/h), where the sum runs over 
all periodic orbits. The amplitude contains information about the stability of a certain 
orbit fj,, with being its action, and its period. All A^, S^, and depend on E and 
X. The coefficient t] acts as a regulator of the infinite summation. Below, our approach will 
cover only the range where the classical perturbation theory is valid, which, nonetheless, can 
correspond to large changes in the spectrum, depending on h |jl8|. With that in mind, we 



assume that a shape distortion will only change the action of a periodic orbit by >S'^(^) ~ 
S^{Xo) + Q^{X — Xq), with = dS^i^Xo) / dX. The change in the orbit itself is assumed to 
be a higher order effect, allowing us to neglect variations in as compared to Qp,/h. After 
some standard manipulations, the diagonal term of can be written as 

Fi-^{X- E) = ^(^Y: MQI {-J-Q^^ - f ^m) ) • (9) 

Two further steps allow us to establish the connection between Eq. (|) and the asymptotic 
part of level velocity correlator. Due to the chaotic nature of the system, the longer the orbits 
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become, the more they explore ergodically the phase space. One can assign a scale Terg to the 
time when periodic orbits with period T > T^rg begin covering uniformly the phase space. 
For T > Terg the Hannay-Ozorio de Almeida sum rule is valid, and the sum in Eq. 
can be enormously simplified |2^, since 0(T) = (j^fi \^tJ.\'^^{T ~ T^)^ — ^ l/27r^|T|. The 



average (. . .) is now taken over a small window 5T which, nevertheless, contains a large 
number of periodic orbits. Notice that this expression specializes our results, since it is 
stated in a form suitable only for time-reversal symmetric systems. (It is trivial to extend 
our results to systems where this symmetry is absent.) In order to proceed analytically, one 
has to extrapolate the Hannay-Ozorio de Almeida sum rule beyond its range of applicability, 
i.e., for T < T^rg [0]. This approximation is not controllable, but its accuracy is supported 
by numerical evidence 0. Hence, 

F^ia9^X;E) ^ fdT<j,{T) ^Qjexp (^Qm^))^^ , (10) 

with (■ ■ ■)sT now indicating an average over periodic orbits within a window 6T, but at a 
fixed energy E. This average can be easily handled following the ansatz proposed in Ref. 



22 , which reads 



P{Q)= I "^Pf-^7^V ^ith {Q^)sT = a{E)T , (11) 

where P{Q) is the probability density of Q. Thus, a is a system-dependent quantity which 
contains information about the long-time dynamics, therefore probing some of its global 
(average) properties. For billiards in the semiclassical regime, the energy dependence of a 
follows straightforwardly from the fact that to a trajectory of length L corresponds an action 
S = {2mEf''^L, yielding 



{Q^iT))sT = 2mE dt' dt 
Jo Jo 



' dL{t') dL{t) \ 
^~dX djT/^^ 



oo 



2mET / dt C{t) . (12) 







This relation is similar to that obtained in Ref. for the variance of the winding number of 
classical trajectories. Here, dL{t) gives the rate in time by which the length of a trajectory 
changes under a variation dX in the billiard shape. For a chaotic system one expects not 
only that {dL(t)/dX)sT = (because different trajectories are uncorrelated) , but also that 
the autocorrelator function of dL{t)/dX decays sufficiently fast to make the integral over 
C{t) finite. The magnitude of / dtC{t) is determined solely by the velocity of the particle 
and the billiard area A. By a proper rescaling one can write ce{E) = {SmE^AY^'^K, where 
K is a classical, dimensionless quantity computed for an unit area billiard for trajectories 
with unit velocity. The energy dependence of a is markedly different from the case where 
the magnetic field plays the role of the external parameter. Notice that a{E) matches the 
largest contribution to C^(0) found in our numerical fit [Eq (^]. 

A this point the semiclassical interpretation of universality becomes apparent: The clas- 
sical long-time dynamics of chaotic systems has some universal properties which are the 
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only necessary elements to evaluate density-density correlators in the semiclassical approx- 
imation. In addition to that, from the classical dynamics one can obtain system-specific 
information which helps to understand relevant scales for different experimental situations. 

After replacing the average of the term involving by an integral over the distribution 
defined in Eq. (|11]), we obtain 

<ri/2j0;i?) = ^, (13) 



Equation (|1^) gives a semiclassical estimate for C„(0). Upon the usual rescaling [Eq. (^)], 
we get 

which, as expected, exhibits the correct large-x asymptotic behavior for time-reversal sym- 
metric systems |]15[, namely, /(x; E) — > —2/ (vr^x^) for x ^ 1. 



The classical evaluation of k makes it possible to predict C^(0). For the conformal 
billiard, for instance, this implies in a simulation of classical trajectories in slightly different 
geometries. We emphasize that in order to give an accurate prediction for k in experiments, it 
is necessary to know very precisely the actual geometry of the quantum dot in consideration. 



V. CONCLUSIONS 

In this paper, we have argued that the universal form of the parametric correlator of 
conductance peak heights indicates the chaotic nature of the electron dynamics in quantum 
dots in the Coulomb blockade regime. In the absence of magnetic fields, the simplest pa- 
rameter to vary experimentally is the shape of the quantum dot. Whereas RMT provides 
the framework for obtaining the universal form of the correlation function, the nonuniver- 
sal correlation scale can be understood in simple semiclassical terms. This scale is rather 
sensitive to the geometry of the dot and the Fermi energy. 

We have compared a theoretical curve obtained from numerical simulations of random 
matrices with the exact correlator obtained for the conformal billiard after averaging over 
energy and shape deformation. The agreement found was good, given the limitations im- 
posed by the size of our data set. We also found that the semiclassical theory can explain 
the energy scaling of C,,(0) on very generic grounds - shedding some light on how to make 
quantitative predictions about nonuniversal characteristics of such systems from a different 
(purely classical) perspective. We remark, however, that shape variations are harder to treat 
classically, since they alter trajectories strongly. To obtain fair estimates of the scaling one 
needs to evaluate {Q'^)st with great numerical accuracy and use very small increments in 
shape deformation. Further investigation in this line is required to make a firm statement. 
It is interesting to note, though, that the energy dependence of the scaling is distinct for 
variations of boundary shape or external magnetic field. 
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FIGURES 




FIG. 1. The rescaled level velocity correlator Cv{x) as a function of shape deformation. The 
full line is the RMT result for the GOE matrix diagonalizations while the points correspond to 
exact numerical calculations for the conformal billiard. 
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FIG. 2. The conductance peak height correlator Cg{x). The symbols follow the same conven- 
tions of Fig. 1. 
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